Time-dependent thermal effects in GRB afterglows 
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Time-dependent thermal effects should accompany standard non-thermal afterglows of GRB when gamma-rays 
pass through inhomogeneous surroundings of the GRB site. Thermal relaxation of an optically thin plasma is 
calculated using time-dependent collisional ionization of the plasma ion species. X-ray emission lines are similar to 
those found in the fading X-ray afterglow of GRB 011211. Thermal relaxation of clouds or shells around the GRB 
site could also contribute to the varying late optical GRB afterglows, such as in GRB021004 and GRB030329. 



1. Time-dependent X-ray emission from a 
thermal plasma 

Determination of physical parameters of as- 
trophysical sources with plasma densities up to 
10 11 — 10 12 cm -3 should include time-dependent 
ionization processes. The cooling time for a 
completely ionized plasma composed of electrons 
and ions with charge Z (due to free-free radia- 
tion alone; allowing for the recombination emis- 
sion decreases this time severalfold) is t c ~ 2 x 
10 15 [s]n eV /T/10 8 /Z In real GRB, the duration of 
gamma-ray emission is 1-100 s and can be much 
smaller than the plasma cooling time. 

X-ray lines in the GRB 011211 afterglow. 
XMM-Newton observations of GRB011211 after- 
glow 11 hrs after the GRB are best fit by a 
thermal plasma continuum with T of several (2- 
4) keV and emission lines (identified with highly 
ionized metal species, except for iron) faded away 
over ~ 10 4 s. In contrast to purely geometri- 
cal interpretation of these observations (which re- 
quires a very high density shell with n e ~ 10 15 
cm), our time-dependent approach identifies this 
time with thermal relaxation time of a plasma 
heated by passing gamma-rays up to a temper- 
ature of several keV (deposited energy eo ~ 10 
keV per nucleon) . This time implies plasma den- 
sity n e ~ 10 11 — 10 12 cm' 3 . Then the observed 
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X-ray emission measure EM = n^V ~ 10 69 cm 3 
means V ~ 10 47 - 10 45 cm 3 . 

Need for clouds. The optical thickness w.r.t. 
Thompson scattering nar < 1 requires a clumpy 
surroundings with the number of clouds N c i ~ 
10 6 of size Id ~ 10 13 cm and the total volume pro- 
viding the required emission measure. (This is an 
upper limit: the realistic cooling function yields 
EM ~ 4 x 10 68 cm~ 3 for normal cosmic abun- 
dance and EM ~ 10 65 cm -3 for plasma with- 
out H and He and other elements in solar ratio). 
The ISM with such cloudy properties is common 
around e.g. late AGB stars Sub-AU clouds 
of the density ~ 10 12 cm~ 3 were also involved by 
PI to explain the observed short time-scale vari- 
ability on GRB light curves in the external shock 
scenario. 

The energy deposit into surroundings. The re- 
quired energy deposit per nucleon eo ~ 10 keV 
limits the external space where such a medium 
can be heated by gamma rays to the distance 

d = - 2 x 10 17 cm where £L = 5 x 10 52 

erg is the total isotropic GRB energy, o 1 k O.Ictt 
is the effective gamma-ray deposition cross sec- 
tion for heating the environment per nucleon 
@|. This distance is independent of the as- 
sumed isotropy/beaming of the GRB emission. 
Again, this distance is strikingly similar to what 
is needed in the external shock scenario [3]. 

GRB jet opening angle 8 should be specified 
to calculate the total volume of the circumburst 
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matter illuminated by GRB. Assuming time de- 
lay t[ = 10 4 s to be entirely due to thermal re- 
laxation and not to geometrical retardation yields 
the relation d(l - cos 6) = (ci/)/(l + z) ~ 10 14 cm 
For small 9, d — 2 x 10 17 cm and z = 2.14 we 
get 9 ~ y/2cti/{\ + z)d ~ 0.05. Such jet colli- 
mation angles 29 ~ 0.1 are common in GRB 5 . 
The total illuminated volume is thus V ~ # 2 d 3 rj 
3 x 10 49 cm -3 which implies the cloud filling fac- 
tor / - 10 45 /10 49 ~ 10~ 4 . So the physically ac- 
ceptable structure of clumpy ISM can actually be 
the source of the observed thermal X-rays from 
GRB 011211. For example, such a medium can 
be produced if several hundred late-type stars 
with strong winds fall within the GRB cone. Un- 
usual in our Galaxy, this could possibly take place 
in star-forming galaxies at large redshifts where 
GRB are observed. 

Modeling the X-ray spectrum. We use time- 
dependent code to compute the state of ionization 
and radiation of thermal plasma which include 
the basic elementary processes (collisional ioniza- 
tion, autoionization, photorecombination, dielec- 
tronic recombination, ion charge exchange) and 
computes free-free, free-bound, bound-bound and 
two-photon emission for ions of all types. The pa- 
rameters include time-dependent temperature T, 
density p and chemical composition. The calcu- 
lated spectral lines are shown in Fig. ^ (see 6 
for more detail). 

2. Time-dependent thermal effects in opti- 
cal GRB afterglows 

The GRB explosion in an inhomogeneous 
medium is a generic picture which is responsi- 
ble for the observed variety of effects in X-ray 
and optical afterglow light curves (see |7j for pio- 
neer study). We consider two possible types of 
surroundings: a cloudy, patchy medium and a 
shell-like structure around the GRB. The second 
case could be identified with an earlier SN ex- 
plosion or intensive stellar wind of the underly- 
ing SN-type event. In both cases, a fraction of 
the GRB energy is deposited into the surround- 
ings. This energy is thermalized and radiated 
away. If plasma of the clouds is optically thin, 
it rapidly cools down. The time-dependent ef- 
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Figure 1. Time evolution of X-ray spectrum (in 
the observer frame at the GRB redshift z — 2.14) 
for different time moments. Upper panel: solar 
composition; Lower panel: metal-rich plasma 
(without H and He). 



fects from the "blue" (approaching the observer) 
gamma-ray beam are expected on the plasma 
cooling timescale t c after the GRB, while the re- 
verberation effects could be observed at different 
wavelengths later on from the "red" (receding) 
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jet on timescales t 2 ~ R(t)/c. Here R(t) is time- 
dependent distance of the photosphere (if formed) 
in the medium from which low-energy photons 
are generated. The maximum time the time- 
dependent effects from the "red" beam will be 
observed is determined by the limiting distance 
d and also by possible hydrodynamical effects. 
So the evolving "SN-like" features in the optical 
spectra observed in the afterglow of GRB 030329 
several days - a month after the burst (also called 
SN 2003dh 8 9 ) could be a manifestation of such 
effects and might have no direct relation to SN 
explosion. 

To illustrate how the cooling plasma forms a 
photosphere, we present the results of LTE calcu- 
lations of plasma opacities as a function of decay- 
ing temperature for a fiducial plasma density p = 
10 -14 g/cm 3 with normal cosmic abundance (Fig. 
EJ. The scattering opacity is ktP ~ 3 x 10~ 15 
cm^ 1 and is independent of temperature until re- 
combination at low temperatures lowers it down 
to ~ 10~ 19 cm -1 . These calculations indicate 
that up to limiting distances d ~ 10 17 cm (as 
appropriate to the case of GRB 030329), such 
a medium remains optically thin for scattering, 
but a photosphere start forming for hard photons 
when temperature drops below ~ fO 6 K. 

To have an idea what kind of optical events can 
be expected from clouds illuminated, we apply 
the multi-group radiation hydrocode STELLA 
[TU| to calculate UBV-light curves (in the GRB 
rest-frame). We assume an energy deposit AE 
thermalizcd in a homogeneous cloud (p = 10 -14 
g/cm 3 ) of mass M over the GRB duration time 
At. No equilibrium between photons and plasma 
is assumed. These calculations are similar to our 
"mini-supernova" model but the difference 
is that here the cloud is transparent for photons 
so the energy is thermalized within the entire vol- 
ume of the cloud. Fig. Olshows the light curve for 
a large (R = 6 X f 15 cm) cloud with a total mass 
of - 5M , AE k 6 x 10 49 erg, and St = fOOs. 
The cloud is heated up to a temperature of 10 4 K. 
The optical emission from the cooling cloud con- 
tinues for more than a month. For comparison, U 
(filled triangles), B (open squares), V (open cir- 
cles) data of a typical plateau type 11 SN f 969L 
are shown. A very bright optical outburst (the 
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Figure 2. Scattering opacity of cooling plasma 
{p = 10~ 14 g/cm 3 ) with normal abundance. 



peak absolute magnitude My ~ —21.5™) is ex- 
pected shortly after the energy deposit (the light 
travel time correction will somewhat smear the 
peak on the time scale ~ R/c.) This effect can 
be compared with the enigmatic outburst seen in 
the R afterglow of GRB 030329 - 1.7 day after 
the GRB [E] . 

A different picture is obtained if we consider a 
smaller homogeneous cloud with R — 1 A x 10 13 
cm and p = 10 -13 g/cm 3 (i.e. like those required 
to explain non-stationary X-ray emission lines in 
GRB 011211). An energy deposit of 10 43 erg into 
the inner AM = 6 x 10~ 7 M Q over At = 10 s heats 
up the cloud to T - 40000 K. This yields (Fig. 

EJ M opt 5 m ... - 7 m . A collection (10 5 - 10 6 ) 

of such clouds (patchy presupernova wind?) at 
d ~ 10 17 cm would naturally produce a reverber- 
ating optical bump with M op ~ — 17™. .. — 19 m on 
the optical light curve of the canonical power-law 
GRB afterglow several weeks after the GRB. It is 
still to be explored if such clouds can mimic spec- 
trum of a SN Ic event spectroscopically identified 
in GRB 030329 [HE]. 
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Figure 3. The "mini-SN" effect from a large ho- 
mogeneous cloud with R = 6 x f 15 cm after en- 
ergy deposit AE w 6 x 10 49 erg. SN 1969L is 
shown for comparison. 



Figure 4. Thermal relaxation of a small ho- 
mogeneous cloud with R = 1.4 x 10 13 cm af- 
ter thcrmalizing AE = 10 43 erg inside the inner 
AM = 6 x 10~ 7 M ra . 



3. Conclusions 

Time-dependent thermal plasma effects from 
GRB surroundings should manifest in X-ray and 
lower energy afterglows of a GRB twofold: 

(1) Shortly after GRB (on the plasma cooling 
time scale determined by density, chemical com- 
position and the energy deposited over the GRB 
emission time) thermal effects may be produced 
by plasma heated up by the "blue" GRB jet. The 
fading X-ray emission lines in GRB 011211 may 
provide an example. 

(2) At later times, variable thermal optical fea- 
tures can be observed superimposed on the under- 
lying non-thermal afterglow continuum due to the 
contribution from cooling clouds / shells heated by 
the receding ("red") GRB jet on timescales lim- 
ited by the most remote clouds d. Optical after- 
glows of GRB 021004 and GRB 030329 deviat- 
ing from pure power-law can provide the exam- 
ple. The effective photosphere moving across the 



heated region can mimic the SN-like effects ob- 
served in some GRB. 
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